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Abstract 

The partial transpose of a block matrix M is the matrix obtained 
by transposing the blocks of M independently. We approach the no- 
tion of partial transpose from a combinatorial point of view. In this 
perspective, we solve some basic enumeration problems concerning 
the partial transpose of permutation matrices. More specifically, we 
count the number of permutations matrices which are equal to their 
partial transpose and the number of permutation matrices whose par- 
tial transpose is still a permutation. We solve these problems also 
when restricted to symmetric permutation matrices only. 
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1 Introduction 



The partial transpose (or, equivalently, partial transposition) is a linear alge- 
braic concept, which can be interpreted as a simple generalization of the usual 
matrix transpose. In the present paper, we consider partial transpose from 
a combinatorial point of view. More specifically, we solve some enumeration 
problems concerning the partial transpose of permutation matrices. 

Even if this notion is a natural one, to the knowledge of the authors, it has 
never been directly studied by the linear algebra community. On the other 
hand, partial transpose is an important tool in the mathematical theory of 
quantum entanglement. For this reason, partial transpose appears often in 
works contextual with quantum information theory. We will spend a few 
paragraphs on this, just for taking a snapshot of the scenario in which this 
notion arises. 

BruB and Macchiavello [3] give an excellent explanation of the meaning of 
partial transpose in quantum information theory. Its primary use is material- 
ized in the so-called PPT- criterion, where "PPT" stands for Positive Partial 
Transpose. The criterion, firstly discovered by Peres [13] and the Horodeckis 
|TOj (see also [12]), is as follows: if the density matrix (or, equivalently, the 
state) of a quantum mechanical system with composite dimension pq is en- 
tangled, with respect to the subsystems of dimension p and g, then its partial 
transpose is positive. The converse of the implication is not necessarily true. 
However, under certain restrictions, for example, when the dimension of the 
density matrix is six, the PPT-criterion is necessary and sufficient. 

There is a number of problems suggested by the PPT-criterion. In partic- 
ular, in order to shed light onto the structure of the set of density matrices, 
it would be important to characterize those for which the criterion is valid. 
An open question of practical importance is to prove or disprove that certain 
states, which are said to be non-distillable, have positive partial transpose. 
However, there is strong evidence that there exist non-distillable states with 
negative partial transpose, which would be then called NPT-bound entan- 
gled states. Regarding this topic, see the important references [SI El, or [1], 
for an account on recent discussions. 

Looking at the notion of partial transpose from the combinatorial point 
of view is an appealing topic, because it has the potential to uncover patters 
in the set of density matrices and indicate connections with other mathemat- 
ical objects, and this may turn out to be helpful in understanding physical 
properties. As a matter of fact there have been a number of recent papers 
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considering entanglement in discrete settings (see, e.g., [Tl [SI fTT]). 

Here we state and solve some basic enumeration problems involving par- 
tial transpose of permutation matrices. Permutations appear in fact to be 
a simple, yet a rich territory to explore. Enumeration is a good first step 
towards the quantitative understanding of the structure of a set. 

In particular, we count the number of permutations matrices which are 
equal to their partial transpose and the number of permutation matrices 
whose partial transpose is still a permutation. We solve these problems also 
when restricted to symmetric permutation matrices only {i.e., induced by 
involutions). 

Apart from considerations related to symmetry, given that symmetry of- 
ten predisposes to relations between different combinatorial objects, a further 
reason to look at involutions comes from [1]. A permutation matrix associ- 
ated to a involution without fixed points can be seen as the adjacency matrix 
of the disjoint union of matchings and self-loops. Since the combinatorial 
Laplacian of any graph is a density matrix after appropriate normalization 
PP, counting the number of involutions whose partial transpose is a permu- 
tation, is equivalent to count the number of these states with positive partial 
transpose. However, the PPT-criterion is not sufficient also for this extremely 
restricted class. There actually are disconnected graphs whose Laplacian is 
entangled even if its partial transpose is positive [9]. 

The organization of the paper is as follows. In the next section, we give 
the required definitions and formally state our problems. Section [3] deals 
with permutations whose partial transpose is a permutation; Section HI with 
permutations equal to their partial transpose; Section [5l with involutions 
whose partial transpose is a permutation 

2 Definitions, statements of the problems and 
examples 

The following is a formal definition of the partial transpose of a matrix: 

Definition 1 Let M he an n ^ n matrix with real entries. Let us assume 
that n = pq, where p and q are chosen arbitrarily. Under this assumption, 
we can look at the matrix M as partitioned into p^ blocks each one qxq. The 
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partial transpose of M, denoted by M p, is the matrix obtained from M, by 
transposing independently each of its p^ blocks. Formally, if 



( B 



1,1 




M = 



\ • • • Bp^p J 



then 



(tot ... »T 




\ • • • 



where Bfj denotes the transpose of the block Bij, for 1 < i, j < p- 



Notice that, by taking the adjoint Bjj, instead of the transpose Bfj, the 
notion of partial transpose can be easily extended to matrices with complex 
entries. This is something which we will not need here. Note that we have 
defined partial transpose with respect to the parameter p. We could have 
also defined partial transpose with respect to the parameter q, by treating 
the blocks of M as the entries of a p x p matrix. Formally, 



That is, the block Bij in M is the block Bj^i in M^", for all 1 < i,j < p. 
The term "partial transpose" also indicates the actual operation required to 
obtain the matrix partial transpose as defined here. 

We will consider partial transpose of permutation matrices. Let us recall 
that a permutation matrix of size n is an n x n matrix, with entries in 
the set {0, 1}, such that each row and each column contains exactly one 
nonzero entry. A permutation of length n is a bijection tt : [n] — > [n] , where 
[n] — {1. 2, ...,n}. Given an n x n permutation matrix P, there is a unique 
permutation tt of length n associated to P, such that 7r(i) = j if and only 
if Pij = 1. Let us denote by Sn the set of all n x n permutation matrices. 
With an innocuous abuse of notation, we write Sn also for the set of all 
permutations of length n. 

In standard linear notation, a permutation vr G S'„ can be written as a 
word of the form 7r(l)7r(2)...7r(n). It may be interesting to point out that a 
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permutation and its partial transpose share a common property, related to 
the sum of the row indices. The cells in the table below contain ordered pairs: 
each left element of the pairs is a permutation n G S^, each right element is 
the ordered list of the row indices of the one entries in the matrix P^^: 



1234, 1234 


1243, 1243 


1324, 1414 


1342, 1432 


1423, 1441 


1432, 1432 


2134,2134 


2143,2143 


2314,4114 


2341,4123 


2413,1414 


2431, 1423 


31 12, 231 1 


31 12, 2323 


321 L 3211 


3211.3223 


3112,3112 


3121.3121 


4123,2341 


4132,2332 


4213,2314 


4231,2323 


4312,4312 


4321,4321 



Example 2 If 



then 
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The matrix P is induced by the permutation tt = 
the row indices of the one entries in P^^ is 2323. 



3124. The ordered list of 



For every permutation n & Sn, where n = pq, we have 

n n 

= n(n +l)/2. 



i=l 



This is straightforward. Let the {ap + i,b{a,i)p + j{a,i))-th entry of P be 
equal to 1. Then b{a, i) runs p times over 0, . . . , g — 1 and j(a, i) runs q times 
over 1, . . . ,p. Thus, 



and 



p+l 
2 
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Therefore 



^ap + j{a,i) =^(2) +^ 



p+1 
2 



n + 1 
2 



which vahdates Eq. ([T]). 

Let us recall that a permutation matrix P is said to be a involution 
if P = P^ and P is not the identity matrix. We will solve the following 
problems: 

Problem 3 Count the number of permutation matrices P G Spq such that 

Example 4 When p = g = 2, we have all together 12 matrices P & S4 such 
that G 5*4 . Among these, 8 are the block-matrices of the forms 



* 
* 



and 



* 
* 



(2) 



The remaining 4 matrices are 



/ 1 








o\ 




/o 


1 





o\ 




f 





1 


\ 




( ° 








1 \ 











1 










1 










1 










1 

















1 
















1 




1 
















1 










1 





0) 












^) 












1 J 









1 


0/ 



(3) 



Problem 5 Count the number of permutation matrices P G Spq such that 



Example 6 When p = g = 2, we have all together 10 matrices P E S4 
such that P^'^ = P. Among these, 8 are the block matrices in Eq. ([2]). The 
remaining 2 matrices are the first and the third matrix in Eq. ([3]). 

Problem 7 Count the number of involutions P G Spq such that P^p G Spq. 
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Example 8 When p = g = 2, we have all together 8 involutions P G 6*4 
such that P^^ = P: 
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3 Permutations whose partial transpose is a 
permutation 

We begin with Problem [31 For a permutation matrix P G Spq, let us de- 
note by Bi^j the block located in the i-th row and j-th column. Let further 
Ajj, Bij C [g] = {1, 2, . . . , g} be the sets of relative row indices and column 
indices of the I's in the block Bi^j. For example, given 

/ 1 \ 

10 

10 ' 
V 1 / 

we have 

= {2},Ai,2 = {1}, ^2,1 = {1},A2,2 = {2}, 

and 

= {1}, fii,2 = {2}, 52,1 = {2}, ^2,2 = {1}. 

Clearly, has the same cardinality as Pjj, which we denote by ri^j. For 
fixed Ai^j and -Bjj, we have r^j! ways to place I's in Bi^j. Therefore, the 
number of required matrices equals the number of Ajj, i^j j's multiplied by 
fl^^. Tj ,,!. At this stage, we impose the required constraints on and Bi j. 
We know that P is a permutation matrix if and only if 

Aij n Ai^k = 0, for every i, j, k with j ^ k, (4) 

Bi j n Pfej = 0, for every k with i k, (5) 



7 



and 

p 

|jA„ = [g], forz = l,2,...,p, (6) 

i=i 
p 

|j5„ = [g], fori = 1,2,..., p. (7) 

i=l 

We need that P^*' is also a permutation matrix. Therefore, we have 

Aij n Akj = 0, for every k with i k, (8) 
n -Bj^fc = 0, for every k with j 7^ fc, (9) 

and 

p 

[jA,, = [q], for J = 1,2,..., p, (10) 

j=i 
p 

|j5„ = [g], for 2 = 1,2,..., p. (11) 

i=i 

Let 

p p 

A^ = ^Ai^^^ and 5^ = P|Si,^^, WireSp, 

i=l i=l 

By Eqs. dH)-©, we know that 

= U A^, S,,, = U 5^. (12) 

From (jl]) and ([S]), we can then write 

At, n Aa- = n Ba- = 0, for every n,a E Sp with i\ ^ a. (13) 
Furthermore, 

U = U = [^]- (14) 

Conversely, given two set partitions {A^} and {B^,} of [q], satisfying Eqs. 
f|T3|) and f|T^ . we may define Aij and by Eq. f|T2l) . One can easily check 
that Eqs. P|)- f|TT|) hold. The only restriction on the A^'s and the -Bjr's is that 
the cardinalities of Aij and Bij should be the equal. Let a.,, and b-„ denote 
the cardinalities of and i?^, respectively. On the basis of the above lines, 
we can state the following result: 
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Theorem 9 Let Z{p,q) the number of permutation matrices P G Spg such 
that P^p G S„a- Then 



pq 

2 



where the sum runs over all aT^^h^^ G Z. 

The following corollary shows a neat expression for the special case P G 

Corollary 10 The number of permutation matrices P G 6*25 such that P^"^ G 
S2q is 

Z(2,g) = g!(g+1)!. 

The pattern avoidance language is now a standard tool for characterizing 
classes of permutations (see [16j). It would be natural to find a character- 
ization of the set of permutations given in Theorem [9] in terms of pattern 
avoidance. 



4 Permutations equal to their partial trans- 
pose 

We now focus on Problem[5l We then ask that Bij = Bfj. Hence, Aij = Bij. 
Additionally, given Aij, the number of ways to put I's in the block Bij is 
exactly the number of involutions of length q, which we denote by I{q)- It is 
well-known that (see, e.g., [15\, Example 5.2.10) 

j even 

and 

/(g + l) = /(g) + g-/(g-l). 

With the same analysis carried on for Theorem [9l we can directly obtain the 
number of desired matrices: 
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Theorem 11 Let Ze{p, q) the number of permutation matrices P G Spq such 
that P = P^p. Then 



where the sum runs over all G Z, with n G Sp. 

When taking p = 2, the number of permutation matrices is given in the 
next corollary: 



Corollary 12 The number of permutation matrices P G 5*2^ such that P = 
P^^ zs 



5 Involutions whose partial transpose is a per- 
mutation 

In this section, we present a solution of Problem [71 Let P be the involution 
defined by the ordered pairs {aq+i, bq+j), where < a, 6 < p—l, 1 < i, j < q 
and (a, i) 7^ {b,j). Note that the partial transpose keeps fixed the I's in the 
diagonal. So, the only possible permutation matrices after partial transpose 
would be the identity matrix Id or P itself. In the first case, we must have 
P = Id, since we get back to the original matrix by applying twice the partial 
transpose operation. Therefore, we only need to consider the second case, 
that is, when P remains invariant under partial transpose. Notice that the 
{aq + i,bq + j)-th and the {bq + j, aq + i)-th entry of the permutation matrix 
are I's. After partial transpose, the {aq+j, bq+i)-th and the {bq+i, aq+j)-th 
entry are I's. Thus we have 




(17) 




r=0 



{aq + i,bq + j) 
{bq + j, aq + i) 



{aq + j, bq + i), 
{bq + i, aq + j). 



or 



{aq + i,bq + j) 
{bq + j, aq + i) 



{bq + i, aq + j), 
{aq + j,bq + i). 



10 



That is, i = j OT a = b. Hence, the desired involutions are of type (ag + 
i, aq + j), with i 7^ j, or, of type {aq + i,bq + i), with a ^ h. We can then 
state the following fact: 

Theorem 13 Let Zt{p,q) be the number of involutions P G Spg such that 
P^'' G Spg, or, equivalently, P^p = P. Then 



The numbers Zt{q + 1, q)/2 are called octagonal pyramidal numbers, and 
count the ways of covering a 2g x 2g lattice with dominoes with exactly 
2 horizontal dominoes (p3], Seq. A002414). The numbers Zt{q,q) count the 
possible rook moves on an a g x g chessboard ([HIH], Seq. A002414). 

To conclude this section, even if these are simple facts, it may be clarifying 
to remark the following: 

Proposition 15 The following statements hold true for all p and q: 

• Z{p,q) = Z{q,p); 

• In general, Ze{p,q) 7^ Ze{p,q); 

• Zt{p,q) = Zt{q,p). 

Proof. While the second point is obvious, the other two can be verified by 
the following bijection. Suppose that the {ap + i, b{a, i)p + j{a, i))-th entry of 
P is 1. Then let the {{i - l)q + (a + 1), {j{a, i) - l)q + (fe(a, i) + l))-th entry 
of P' be 1. If the partial transpose of P is a permutation, then ap + j{a,i) 
and b{a, i)p + i run from 1 to n, for < a < g — 1, 1 < z < p. Thus, 
{i — l)g + (6(a, i) + 1) and j(a, i) — l)g + (a + 1) run from 1 to n also. This 
implies that the partial transpose of P' is a permutation. ■ 
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Corollary 14 The following statements hold true: 



• Zt(g + l,g)=g(g + l)(2g-l); 
. Z,(g,g) = 2(g3-g2). 
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